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ABSTRACT 

A dense ionized cloud of gas has been recently discovered to be moving directly 
toward the supermassive black hole, Sgr A*, at the Galactic Center. In June 2013, 
at the pericenter of its highly eccentric orbit, the cloud will be approximately 3100 
Schwarzschild radii from the black hole and will move supersonically through the am- 
bient hot gas with a velocity of v p ~ 5400 km s _1 . A bow shock is likely to form in 
front of the cloud and could accelerate electrons to relativistic energies. We estimate via 
particle-in-cell simulations the energy distribution of the accelerated electrons and show 
that the non-thermal synchrotron emission from these electrons might exceed the qui- 
escent radio emission from Sgr A* by a factor of several. The enhanced radio emission 
should be detectable at GHz and higher frequencies around the time of pericentric pas- 
sage and in the following months. The bow shock emission is expected to be displaced 
from the quiescent radio emission of Sgr A* by ~ 33 mas. Interferometric observations 
could resolve potential changes in the radio image of Sgr A* at wavelengths < 6 cm. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active - 
Galaxy: center 

1. Introduction 

Recent sub-millimeter observations revealed a dense ionized cloud of gas known as G2, rapidly 
approaching Sgr A*, the black hole at the Galactic Center (Gillessen et al. 2012). The cloud is 
on a highly eccentric trajectory, with a 2011 distance from the black hole of 1.8 x 10 16 cm. The 
pericentric passage, which is expected to occur in mid 2013, will bring the cloud within 4 x 10 15 cm 
from the supermassive black hole. Given the mass of the black hole, which is determined through 
observations of nearby stellar orbits to be M = 4.3 x 10 6 M@ (Ghez et al. 2008; Gillessen et al. 2009), 
this pericentric distance is only R p = 3100 Rs, where the Schwarzschild radius Rs = 1.27 x 10 12 cm. 

The accretion flow around the black hole extends to the Bondi radius ~ 10 5 i?s (e.g., Yuan, 
Quataert & Narayan 2003) and powers the multiwavelength emission observed from it. The flux at 
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1 GHz is ~ 0.5 Jy, rising to ~ 4 Jy at 500 GHz, before rapidly declining at higher frequencies. The 
radio emission has been successfully modeled as synchrotron radiation from relativistic electrons, 
either in a radiatively inefficient accretion flow (ADAF) (Narayan, Yi, & Mahadevan 1995; Ozel, 
Psaltis & Narayan 2000) or in a jet (Falcke & Markoff 2000). At the lowest end of the spectrum — 
v ~ 1 — 10 GHz — the radio emission shows flux variability on a time scale of months to years with 
a root mean square amplitude ~ 10%, (Zhao et al. 1989; Falcke 1999; Macquart Sz Bower 2006). 

At its pericentric passage, the gas cloud will interact with the accretion flow, and this may 
significantly change its dynamics. Here, we show that a bow shock is likely to develop as the 
cloud plows through the hot, tenuous plasma at R p . In §2 we calculate the properties of the bow 
shock and in §3 we estimate the energy distribution of electrons accelerated at this shock. In §4 
we calculate the extra radio emission that will result from these accelerated electrons. For likely 
electron energy distributions, the additional emission is ~ 10 Jy at frequencies ~ 1 — 10 GHz. This 
is well above the quiescent emission from Sgr A* and should be easily detectable. We also show 
that the flux increase will be accompanied by significant changes in the spectral index. In §5 we 
summarize our findings and argue that interferometric observations could resolve potential changes 
in the radio image of Sgr A* caused by the interaction of the cloud with the accretion flow. 



2. The Bow Shock around G2 at Pericenter 

We first calculate the dynamics of the interaction of the cloud G2 with the accretion flow 
around the black hole. Since the orbit of G2 is highly eccentric, its velocity at pericenter will be 



v p 



[2GM\ 



1/2 



J = 5400 kms~ 1 . (1) 



The properties of the ambient gas at the pericentric distance R p can be obtained from ADAF 
models of Sgr A*. Within that context and using Yuan et al. (2003) and Xu et al. (2006) as a 
guide, Gillessen et al. (2012) estimated the gas density and temperature at R = R p to be 

n^ 930 ^ L4 X ^° 4 Rs ^j cm~ 3 = 4200 cm~ 3 , T « 10 9 K. (2) 

To estimate the sound speed, we use the fact that the gas in the ADAF is likely to have its Bernoulli 
parameter close to zero (Narayan & Yi 1994). Thus, 



GM 1 9 1 
Be = -— + -v% + -^ + w^0, (3) 



where w = Tp/(T — l)p is the enthalpy per unit mass, T = 5/3 is the adiabatic index, p is the 
pressure, and p the density of the gas. Because the radial velocity vr is small compared to the 



1 Some models suggest a temperature closer to 10 8 ' 5 K, but this small uncertainty is not important for what follows. 
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azimuthal velocity v$, we ignore it. We also assume that v\ is roughly half the Keplerian value, 
which is appropriate for an ADAF. Under these assumptions, we estimate the adiabatic sound 
speed of the gas, c a d = y^Tp/p, to be 

Cad ~ v p /2 ps 2700 kins" 1 . (4) 

The gas pressure is then given by 

P = fPcld w 3.1 x 10~ 4 erg cm- 3 . (5) 

Finally, if we assume that the magnetic pressure is about 10% of the gas pressure, as is typical in 
an ADAF, then we estimate the magnetic field strength in the ambient gas to be ~ 0.03 G. 

There is currently no information on the relative orientation of the orbits of G2 and the ambient 
gas in the accretion flow. However, because the velocity v p of the cloud at pericenter is substantially 
larger than the local gas velocity v<p, we expect the relative velocity at R p to be dominated by v p . 
Thus, G2 will move with a Mach number A4 ~ 2 near pericenter. Using Rankine-Hugoniot jump 
conditions for a non-relativistic shock, we estimate for the shocked gas, 

n sho ck (r + l)M 2 



n (r-l)A4 2 + 2 

Pshock (r + 1) + 2r(M 2 - 1 



2.3, (6) 
4.8. (7) 



V (T + l) 

For a Mach number A4 ~ 2, the magnetic field strength will be roughly doubled by compression 
in the shock, i.e., -B s hock ~ 0.06 G. Similarly, the post-shock temperature will be a factor of a few 
larger than the temperature of the pre-shock gas. Thus, the mean thermal energy of the shocked 
electrons should be close to m e c 2 , i.e., the electrons will be quasi-relativistic. 

Gillessen et al. (2012) estimated the equivalent spherical size of the cloud in 2011 to be 15 mas, 
which corresponds to a physical radius of 1.9 x 10 15 cm. As the cloud approaches the pericenter 
along its highly eccentric orbit, it is expected to be tidally stretched (Gillessen et al. 2012; Burkert 
et al. 2012), primarily along the direction of motion. We estimate the cross-sectional area A of 
the bow shock using the frontal size of the cloud, which will remain approximately constant at 
7r(10 15 cm) 2 . 

Figure Q] shows the Mach number of G2's motion through the ambient medium (assuming that 
the medium is at rest) and the temperature of the medium as a function of time. A bow shock 
will form when A4 exceeds unity, which is expected to last for about a year around pericenter. 
The temperature of the medium, which has a significant effect on particle acceleration (see §3), is 
also strongly peaked around pericenter. Taking the relative velocity between G2 and the ambient 
medium to be equal to v p , and conservatively taking the duration of the bow shock to be 6 months, 
we estimate the total number of shocked electrons to be 

^shock ~Av p t p n = 1.1 x 10 50 . (8) 
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Fig. 1. — Mach number Ai of G2's motion through the hot accretion flow around Sgr A* and 
temperature T of this medium as a function of time. Pericentric passage occurs in year 2013.5. Both 
M. and T peak sharply around pericenter, and a bow shock should be present for approximately 6 
months on either side of pericenter. Since the orientation of the angular momentum vector of the 
accreting gas is not known, for simplicity the gas has been assumed to be at rest. 



3. Acceleration of Electrons in the Bow Shock 

A variety of astrophysical evidence suggests that particles are accelerated efficiently via the 
Fermi and shock-drift acceleration mechanisms in shocks (e.g., Blandford & Eichler 1987). Typ- 
ically, these processes give rise to a non-thermal power-law tail in the energy distribution of the 
particles. The parameters of the bow shock of G2 correspond to an interesting regime that has not 
been well studied. The shock is non-relativistic, but the upstream electrons are quasi-relativistic 
(kT e < m e c 2 ). In addition, the shock has a modest Mach number A4 ~ 2, and the upstream gas is 
fairly strongly magnetized, corresponding to an Alfvenic Mach number Ma ~ 8. 

We studied the acceleration of electrons in the bow shock of G2 by means of two-dimensional 
first-principles numerical simulations, with the particle-in-cell (PIC) code TRISTAN-MP (Spitkovsky 
2005). The simulation setup parallels very closely the one employed by Riquelme & Spitkovsky 
(2011), with the magnetic field lying initially in the simulation plane, oriented at an oblique angle 
with respect to the flow velocity. For computational convenience, we chose a reduced mass ratio 
m p /m e = 100, but we tested that our results remain the same for larger mass ratios, when all the 
physical quantities are scaled appropriately. Specifically, we ran simulations spanning the range 
m p /m e = 25 — 400, fixing the electron temperature (equal to the proton temperature) and the 
shock sonic and Alfvenic Mach numbers, and we measured the time in units of the inverse proton 
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Fig. 2. — Electron energy spectrum just upstream of the shock, normalized to the pre-shock electron 
density. (Left) The temporal evolution of the spectrum, for T e = 10 9 K, M = 2 and Ma = 8. 
With time, the non-thermal component approaches a power-law tail with index p = 2.2 (dashed 
line) and its cutoff energy steadily increases (insert). (Right) Electron energy spectrum at uj c it = 14 
for different upstream temperatures (therefore, different M), with fixed Ma = 8 and fixed flow 
velocity. The regime 10 85 K < T e < 10 9 K is relevant for the accretion flow around Sgr A* at a 
distance R p from the black hole. 

cyclotron frequency u^ 1 . We also checked the convergence of our results with respect to the spatial 
resolution and the number of computational particles per cell. 

We find that at the shock, a fraction of the incoming electrons are reflected backward by the 
shock-compressed magnetic field (Matsukiyo et al. 2011) or by scattering off of electron whistler 
waves excited in the shock transition layer (Riquelme & Spitkovsky 2011). For quasi-relativistic 
electron temperatures, the reflected electrons are fast enough to remain ahead of the shock, resisting 
advection downstream by the oblique pre-shock field. While the electrons gyrate around the shock, 
they are energized by shock-drift acceleration (e.g., Begelman & Kirk 1990) and form a local non- 
thermal population, just upstream of the shock. 

The temporal evolution of this population is shown in the left panel of Figure for parameters 
relevant to the bow shock of G2: T p = T e = 10 9 K ~ 0.2 m e c 2 /k, M = 2 and Ma = 8. At late times, 
the shock-accelerated electrons populate a power-law tail with a slope of p = 2.2 beyond an electron 
Lorentz factor 7 e ~ 2, containing roughly 5% of the incoming electrons. The upper energy cutoff 
of the electron spectrum steadily increases with time, as shown in the left panel insert of Figure 
suggesting that the distribution will asymptote at late times to a power-law with p < 2.2 extending 
to very large values of j e . The counter streaming between the incoming flow and the shock-reflected 
electrons triggers the Weibel filamentation instability ahead of the shock (Weibel 1959, Medvedev 
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& Loeb 1999), with the wavevector perpendicular to the pre-shock fieldo By scattering off of the 
magnetic field generated by the Weibel instability, the shock-reflected electrons are deflected back 
toward the shock, participating in a Fermi-like acceleration process. In the downstream region, 
they populate a power-law tail of similar normalization and slope as the pre-shock spectrum shown 
in the left panel of Figure O 

The distribution of non-thermal electrons is sensitive to the electron temperature ahead of the 
shock. As shown in the right panel of Figure [2j the normalization of the power-law tail is reduced 
by almost one order of magnitude when the electron temperature decreases from T e = 10 9 -ftT (red 
curve) down to T e = W 8 K (blue curve), with the flow velocity staying fixed. If the upstream 
plasma is colder, fewer electrons are reflected back at the shock (Matsukiyo et al. 2011), so a 
smaller fraction of the incoming electrons are injected in the shock-drift acceleration process. In 
the limit of cold upstream plasmas studied by Riquelme & Spitkovsky (2011), electrons are not 
efficiently reflected back from the shock (black curve in Figure [21 for T e = 5 x 10 7 K). In this case, 
the process of Weibel-mediated acceleration described above does not operate and the resulting 
downstream non-thermal tail becomes steeper, with p > 3 (Riquelme & Spitkovsky 2011). 

The acceleration efficiency of 5% and the power-law index of 2.2 that we find for the parameters 
of the bow shock of G2, combined with our estimate for the total number of shocked electrons given 
in Equation (|8|), allows us to write the electron energy distribution as 



In the next section, we calculate the properties of the synchrotron emission that arises from this 
electron distribution. 



The peak of the synchrotron emission from an electron with a Lorentz factor 7 e occurs at a 
frequency 



where we have scaled the result to the expected field strength of 0.06 G in the shocked medium. 
Conversely, we can invert the above relation to infer the Lorentz factor of the electrons that con- 
tribute predominantly to the emission at a particular frequency: 




7e > 2. 



(9) 



4. Expectations for Non-thermal Synchrotron Emission 




(10) 




(11) 



2 We point out that the Weibel mode can only be captured by means of multi-dimensional simulations and was, 
therefore, absent in the one-dimensional experiments of Matsukiyo et al. (2011). 
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The synchrotron power emitted by such electrons is 

iW « 3.8 x 10- (^) \? ergs"' « 1.6 x 10"" (-f^) , GH , ergs"'. (12) 

The synchrotron cooling time is then 

l e m e c 2 (BY 2 -, ( B \~ 3/2 _ 1/2 , 

*~* = W 6800 U06G ) ^ yr ^ 105 (o06G ) yr- (13) 

The cooling time is substantially longer than the duration of the encounter for any value of the 
Lorentz factor j e or, equivalently, any frequency i^GUz of interest for radio or submm observations. 
Thus, all the shocked electrons will contribute to the observed synchrotron emission. 

To estimate the expected spectral flux from the shocked electrons, we use the electron Lorentz 
factor distribution given in Equation © and assume a distance D = 8.3 kpc to the Galactic 
Center. Then, using the standard formula for synchrotron emission from a power-law distribution 
of electrons (Rybicki & Lightman 1979), we obtain 

F "" v> {moT'^^ p = 2 ' 2 - (14) 

This estimate is valid over the range of frequencies at which the synchrotron emission is optically 
thin. The predicted flux is fairly large and should be easily detected over a wide range of radio fre- 
quencies, provided particle acceleration in the bow shock is as efficient as the numerical simulations 
described in §3 indicate. 

At low frequencies, the synchrotron emission will be self-absorbed. Using standard results 
(Rybicki & Lightman 1979), we estimate the source function of the shocked electrons to be 

S u = 3.7 x 1(T 8 ( -j^J 4hz ergcm-V^z-W- 1 . (15) 

Assuming a circular source of radius 10 15 cm at distance D, the limiting synchrotron flux due to 
self-absorption is 

/ B \-l/2 

JW « is [-^s) 4L Jy- (is) 

The quiescent emission from Sgr A* is well below this level at frequencies above a GHz, so self- 
absorption should not interfere with our ability to observe the additional emission. 

In Figure [3j we show the predicted radio emission from shock-accelerated electrons in the bow 
shock of G2. In order to account for potential uncertainties in our estimates of shock parameters 
(§2) and our simulations of particle acceleration (§3), we explore the dependence of the expected 
flux enhancement on (left panel) the number of accelerated electrons and (right panel) the power- 
law index of the electron energy distribution. We also show in Figure [3] the measured quiescent 
flux at different frequencies, with the error bars indicating the degree of variability among different 
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Fig. 3. — Radio emission expected from non-thermal electrons accelerated in the bow shock of G2, 
plotted along with the quiescent emission observed from Sgr A*. The data points are compiled 
from Falcke et al. (1998), Zhao et al. (2003), and Marrone et al. (2008). The solid line in the left 
panel shows the predicted spectrum for our fiducial model with power-law index p = 2.2 (eq. [9J. 
The hatched area bounded by dashed lines corresponds to a factor of three uncertainty each way 
in the number of accelerated electrons, i.e., the acceleration efficiency is varied from 1.7% to 15%. 
The hatched region in the right panel shows the predicted radio flux for power-law indices in the 
range p = 2 — 2.5, keeping the fraction of accelerated electrons fixed at 5%. In both panels, the 
turnover of the flux at low frequencies is caused by synchrotron self-absorption. 



observations. At 1.4 GHz, the quiescent radio flux of Sgr A* is 0.5 Jy, whereas we estimate that 
the flux enhancement could be as large as 10 Jy. The additional synchrotron emission ought to 
be easily detectable at GHz frequencies. Note also that the spectral index is predicted to change 
substantially. 



5. Discussion 



The passage of the recently discovered cloud of gas G2 near Sgr A* presents a unique opportu- 
nity to study the dynamics and properties of hot gas in the vicinity of the black hole at the Galactic 
Center. In this Letter, we showed that a bow shock may form during the pericentric passage of the 
cloud, and we investigated the flux enhancement in the 1 — 100 GHz frequency range that will arise 
as a result of particle acceleration in the shock front. We ran first-principles PIC simulations for 
shock parameters appropriate to the bow shock and thereby obtained realistic estimates of the en- 
ergy distribution of accelerated electrons. Using these results, we calculated the likely synchrotron 
emission from the bow shock and found that the additional flux might exceed the quiescent emission 
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from Sgr A* by up to an order of magnitude at GHz frequencies. This suggests that there is a good 
chance of detecting enhanced radio emission as G2 plows through the ambient hot medium around 
the time of pericentric passage. Since the cooling time of the accelerated electrons is estimated to 
be long, the enhanced emission should continue well after the encounter. 

There are order unity uncertainties in the parameters we have assumed for the bow shock, 
and hence the predictions made in this Letter are not likely to be quantitatively accurate. We 
have allowed for some of these uncertainties while computing the hatched regions shown in the 
two panels in Figure [3l An additional uncertainty is whether or not a bow shock will form in the 
first place. Gillessen et al. (2012) discuss a compression shock moving into the cloud, which will 
inevitably be accompanied by an external shock moving into the ambient medium, the bow shock 
in our model. In most models of G2 (Burkert et al. 2012; Miralda-Escude 2012; Schartmann et al. 
2012; Murray-Clay & Loeb 2012), the cloud retains some level of integrity during its pericentric 
passage and thus is likely to develop an external shock. However, if the cloud is completely shredded 
by Kelvin-Helmholtz or other instabilities before a bow shock forms, then our synchrotron emission 
estimates will no longer be valid. 

Given the pericentric distance of 3100 Rs, which corresponds to a projected angle of ~ 33 mas, 
the bow shock emission should be displaced from the quiescent radio emission of Sgr A* by the 
same amount. At wavelengths < 6 cm, this angular distance is larger than the size of the scattering 
ellipse of Sgr A* (Bower et al. 2006) and ought to be resolved by interferometric observations. The 
scatter-broadening of the radio image of Sgr A* is believed to be caused by a compact foreground 
interstellar cloud that is at least ~ 100 pc from the black hole (Frail et al. 1994). Thus, the 
broadening is unlikely to be affected by any gas stripped from G2 during its pericentric encounter. 

In addition to the bow shock and the associated prompt radio synchrotron emission considered 
in this Letter, it is expected that the cloud G2 will also shed mass as it interacts with the ambient 
hot gas. A likely early signature of the increase in the gas density at a few thousand Rg is a 
change in the observed Faraday rotation above a GHz, which may provide the first estimates of the 
increase in the mass accretion rate. As it moves inwards, this gas will cause the mass accretion rate 
on to the central black hole to be enhanced over a period of many years. Such an increase could 
cause a secular change in the radio flux of Sgr A* on a time scale of ten years to several decades, 
accompanied by changes in the "silhouette" of the black hole that could be monitored by future 
interferometers (Moscibrodzka et al. 2012). 
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